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ABSTRACT 

Millimeter wave (mmWave) systems will likely employ di¬ 
rectional beamforming with large antenna arrays at both the 
transmitters and receivers. Acquiring channel knowledge to 
design these beamformers, however, is challenging due to the 
large antenna arrays and small signal-to-noise ratio before 
beamforming. In this paper, we propose and evaluate a down¬ 
link system operation for multi-user mmWave systems based 
on compressed sensing channel estimation and conjugate 
analog beamforming. Adopting the achievable sum-rate as a 
performance metric, we show how many compressed sensing 
measurements are needed to approach the perfect channel 
knowledge performance. The results illustrate that the pro¬ 
posed algorithm requires an order of magnitude less training 
overhead compared with traditional lower-frequency solu¬ 
tions, while employing mmWave-suitable hardware. They 
also show that the number of measurements need to be opti¬ 
mized to handle the trade-off between the channel estimate 
quality and the training overhead. 

Index Terms — Millimeter wave communication, com¬ 
pressed sensing, achievable rates. 

1. INTRODUCTION 

Millimeter wave (mmWave) communication is a promising 
technology for future cellular systems 00 - Directional pre¬ 
coding with large antenna arrays appears to be inevitable to 
support longer outdoor links and to provide sufficient received 
signal power. The design of precoding matrices, though, is 
usually based on complete channel state information, which 
is difficult to achieve in mmWave due to (i) the large num¬ 
ber of antennas, (ii) the small signal-to-noise ratio (SNR) be¬ 
fore beamforming, and (iii) the different hardware constraints 
which impacts the training signal design 0 - Therefore, de¬ 
veloping low-complexity mmWave channel estimation tech¬ 
niques is crucial for the mmWave system operation. 

To overcome the hardware limitations, analog beamform¬ 
ing solutions were proposed in which rely on control- 
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ling the phase of the signal transmitted by each antenna via a 
network of analog phase shifters. To avoid the need for ex¬ 
plicit channel knowledge, the algorithms in |6|jT0) depend 
on beam training which iteratively designs the analog beam¬ 
forming coefficients at the transmitter and receiver. These so¬ 
lutions, however, have two main disadvantages: (i) they sup¬ 
port only single-stream transmissions (2) their training over¬ 
head scales linearly with the number of users. To support 
multi-stream transmission, |5| proposes an efficient mmWave 
channel estimation algorithm by leveraging the sparse nature 
of the channel and adaptive compressed sensing tools. This 
solution, however, did not solve the overhead scaling issue. 

In this paper, we propose a simple downlink system oper¬ 
ation for multi-user mmWave systems based on compressed 
sensing channel estimation and conjugate analog beamform¬ 
ing. Contrary to prior work in 0®- the training overhead 
of the proposed channel estimation solution does not scale 
with the number of users. Hence, it is of special interest 
for multi-user mmWave systems. Under certain assumptions, 
we characterize a lower bound on the achievable rate of the 
proposed system operation as a function of the compressed 
sensing measurements. Simulation results show that the pro¬ 
posed algorithm requires much less training overhead com¬ 
pared with traditional lower-frequency solutions, while em¬ 
ploying mmWave-suitable hardware. 

We use the following notation: A is a matrix, a is a vector, 
and a is a scalar. |j A||i? is the Frobenius norm of A, whereas 
A 1 , A*, A -1 , are its transpose, Hermitian, and inverse, re¬ 
spectively. E [■] denotes expectation and \ is an indicator. 

2. SYSTEM MODEL 

Consider a mmWave system with a base station (BS) hav¬ 
ing Nbs antennas and (Vrf RF chains as shown in Fig. [T] 
The BS is assumed to communicate with U mobile stations 
(MS’s), and each MS is equipped with Nms antennas. We 
focus on the multi-user beamforming case in which the BS 
communicates with every MS via only one stream. Further, 
we assume that the maximum number of users that can be si¬ 
multaneously served by the BS equals the number of BS RF 
chains, i.e., U < Aj{p. This is motivated by the spatial mul- 





Fig. 1. A BS with RF beamformers and jVrf RF chains com¬ 
municating with the wth MS that employs RF combining. 


tiplexing gain of the described multi-user precoding system, 
which is limited by min (A r RF . U) for AAs > Arf • For sim¬ 
plicity, we will also assume that the BS will use U out of the 
Arf available RF chains to serve the U users. 

On the downlink, the BS applies an A As x U RF precoder, 
F = [fi, f2,.... f/y]. The sampled transmitted signal is there¬ 
fore x = Fs, where s = [si, s 2 , S[/] T is the U x 1 vector of 

transmitted symbols, such that E [ss*] = u~ and P x is the 
average total transmitted power. Since F is implemented us¬ 
ing quantized analog phase shifters, [F] = — 


p]<t>n 


is a quantized angle, and the factor of 


is 


where < 
for power normalization. 

For simplicity, we adopt a narrowband block-fading chan¬ 
nel model 0HE3’ by which the 71 th MS receives the signal 


u 

= H u ^f r 

r=1 


n 7/ 


(i) 


where H„ is the AAis x AAs matrix that represents the 
mmWave channel between the BS and the 71th MS, and 
n u ~ A/"(0, cr 2 1) is a Gaussian noise vector. 

At the uth MS, the RF combiner w„ is used to process the 
received signal r u to produce the scalar 

u 

y u = w*H u ^f r s r + w*n u , (2) 

r =1 

MmWave channels are expected to have limited scatter¬ 
ing ©■ Therefore, and to simplify the analysis, we will as¬ 
sume a single-path geometric channel model [9][T3j. Under 
this model, the channel H„ can be expressed as 

H„ = ATBsAr M sa«aMS (#«) ^bs , (3) 


and each MS are assumed to know the geometry of their an¬ 
tenna arrays. While the results and insights developed in the 
paper can be generalized to arbitrary antenna arrays, we will 
assume uniform arrays in the simulations of Section [5] 

3. PROPOSED DOWNLINK SYSTEM OPERATION 

The proposed downlink operation for multi-user mmWave 
systems consists of two phases: (i) compressed sensing based 
downlink channel estimation and fii) conjugate analog beam¬ 
forming/combining. For the downlink channel training, ran¬ 
dom beamforming and projections are used to efficiently 
estimate the mmWave channel with relatively low training 
overhead thanks to the sparse nature of the channel. One 
main advantage of this technique is that all the MS’s can si¬ 
multaneously estimate their channels. Therefore, the training 
overhead does not scale with the number of users. This is 
contrary to the adaptive channel estimation and beamforming 
design techniques in [5j[6j[9j, which are user-specific. The 
estimated channels are then used to build the analog beam- 
formers and combiners. Extensions to hybrid analog/digital 
precoders are also possible G3 , but our focus in this paper is 
on the evaluation of compressed sensing channel estimation. 

3.1. Compressed Sensing Based Channel Estimation 

Given the geometric mmWave channel model in Q, estimat¬ 
ing the channel is equivalent to estimating the different pa¬ 
rameters of the channel path; namely its AoA, AoD, and the 
complex gain. In this section, we exploit this poor scattering 
nature of the mmWave channel, and formulate the channel es¬ 
timation problem as a sparse problem. We then briefly show 
how compressed sensing can be used to estimate the channel. 

A sparse formulation: Consider the system and mmWave 
channel models described in Section [2] If the BS uses a train¬ 
ing beamforming vector p m , and the uth MS employs a 
training combining vector q„ to combine the received signal, 
the resulting signal can be written as 

yn,m — Bn T H n,m: (4) 

where s m is the training symbol on the beamforming vec¬ 
tor p m , and we use s m = y/P, with P the average power 
used per transmission in the training phase. If the BS em¬ 
ploys Mbs su ch beamforming vectors p m ,m = 1, ...,Mbs, 
at Mbs successive time slots, and the MS uses Mms mea¬ 
surement vectors q n ,n = 1,2, ...,Mms at Mms successive 
instants to detect the signal transmitted over each of the beam¬ 
forming vectors, the resulting received matrix will be |51 


where a u is the complex path gain, including the path-loss, 
with E [|a M | 2 ] = a. The variables 9 U , and cf) u £ [0, 27t] 
are the angles of arrival and departure (AoA/AoD) respec¬ 
tively. Finally, ass (<j> u ) and aMS (8 U ) are the antenna array 
response vectors of the BS and 7ith MS respectively. The BS 


Yms = VPQ H U u P + N, (5) 

where Q = [qi,q 2 , ...,qM MS ] is th e Ams x M m s measure¬ 
ment matrix, P = [pi, p 2 ,..., Pa/ bs ] is the BS N BS x M B s 
beamforming matrix, and N is an Mms x Mbs noise matrix. 































To exploit the sparse nature of the channel, we first vec¬ 
torize the resultant matrix Yms as in {5]| to get 

Yms = VP (P T <g> Q ff ) (agg ((j) u ) ® a M s (0J) + v, 

( 6 ) 

To complete the problem formulation, we assume that the 
AoDs, and AoAs are taken from a grid of Gbs and Gms 
points, respectively. By neglecting the grid quantization er¬ 
ror, we can approximate yms in <[h]) as 

yMS = VP (P T <8> Q h ) ^Agg <g> Ams) Zti + v, (7) 

where the Abs x Gbs matrix Abs and Nms x Gms matrix 
Ams are the dictionary matrices that consist of the column 
vectors ass (</>&), and aMS (^), respectively, with and 61 
the A;th, and £th points of the angle grids. z u is a Gbs Gms x 
1 vector which carries the path gains of the corresponding 
quantized directions. Hence, z u is a sparse vector with only 
1 non-zero element. Note that detecting the column of Abs 
and Ams that corresponds to this non-zero element directly 
implies the detection of the AoA and AoD of the channel path. 

Compressed sensing measurements: Thanks to the 
sparse formulation of the mmWave channel estimation prob¬ 
lem in 0. compressed sensing tools can be leveraged to 
design efficient training beamforming/combining matrices 
1 5jl5 1. Considering the measurement matrix <P = P 7 ® Q H , 
and the dictionary \P = A ISS 0 A ms, one interesting research 
direction is to study the conditions on <P, \P under which the 
support of the sparse vector z u can be recovered with high 
probability and with low training overhead. Leaving this ob¬ 
jective for future work, we will try in this paper to get some 
insights into a sufficient (not necessarily the minimum) num¬ 
ber of measurements, and the relation between this number 
and the achievable rate of mmWave systems. For that, we 
consider the following measurement matrix. 

The BS will design its training beamforming matrix 
P, such that [P] = e j< ™> n where is randomly 

and uniformly selected from the set of quantized angles 

{0, -Jgg,..., Q M tss }• Each MS similarly designs its 

training combining matrix Q, with N^ s angle quantization 
bits. Note that for this design, each entry of the measurement 
matrix $ will also be equal to e^, with the angle ( selected 
randomly from a certain quantized angle set. 

A 0 A/A 0 D estimation: To estimate its channel AoA/AoD, 
each MS u needs to recover the support of its sparse vector 
z u . As the measurement and dictionary matrices <F. \P are 
known for all MS’s, each MS can use sparse recovery algo¬ 
rithms (such as LASSO (T§, Orthogonal Matching Pursuit 
(OMP) [F7j, etc.) to estimate its channel A 0 A/A 0 D. In the 
simulations of Section[5] we adopt OMP for low-complexity. 
In this case, the support of z u , supp(z u ), is determined by 
solving 

supp(z„) = arg max ’P " T> "y MS , (8) 

which directly determines the estimated AoA/AoD ,9 U , 0 U . 


3.2. Conjugate Analog Beamforming 


Given the estimated AoA, each MS u will design its analog 
combining vector such that w„ = aMS • Each MS u 
will also feed the index of its estimated AoD back to the BS 
which needs log 2 Gbs bits. Finally, the BS designs its analog 
beamforming matrix F to match the effective channels (in¬ 
cluding the effect of the combining vectors), i.e., the BS sets 




4. ACHIEVABLE AND EFFECTIVE RATES 


In this section, we evaluate the achievable rate of the proposed 
downlink mmWave system operation in a special case and try 
to get some insights into the relation between the performance 
of mmWave systems and compressed sensing measurements 
in more general cases. 

Consider the system model in Section[2]and the proposed 
conjugate analog beamforming/combining in Section [3T2] For 
tractability, we make the following assumptions 

Assumption 1 All path gains are constants. This is relevant 
to mmWave LOS paths, which are dominant in dense mmWave 
networks m 

Assumption 2 The BS and MS’s employ uniform arrays. 


Assumption 3 The sizes of the grids in 0 are Gbs = Abs. 
Gms = Ams. an d the angle grid points (j)k,9g, of the dictio¬ 
nary \P are the virtual directions satisfying sin = 

7&and** sin(^) = ^M 


Under the assumptions [2|3] the matrices Abs and Ams 
become DFT matrices GD- Further, using the virtual chan¬ 
nel model transformation in (T9), we note that the steering 
vectors of both the actual and estimated AoAs/AoDs (of the 
dictionary tP) are columns of these matrices. Therefore, the 
rate of user u, R u , can be written as 


log 2 1 + 


1 {e u =e„,0„=^} 


Er=i 1 


1 r} ^ WN BS N MS \a u \' 2 


(9) 


Denoting R u = log 2 (l + ^ A B sA M s |a u | 2 ) as the 
single-user rate (without interference), the average achievable 
rate of user u, R u , can be then lower bounded as 


R„ > E 




( a ) p p p 


> Ru 1 — 


U 

Abs 


{B U =6 U , <t> u =$ u }’ 


( 10 ) 

( 11 ) 

( 12 ) 


where (a) follows from the independence between the estima¬ 
tion success event { 9 U = 9 U ,4> U = (f> u } and the single-user 
















Fig. 2. Achievable rates using the proposed system operation Fig. 3. Achievable rates using the proposed system operation 

for different numbers of compressed sensing measurements for different values of channel fading coherence. 


rate given assumption|Tj Now, we note that =g ^ j 
is the probability of correct support recovery of the sparse 
vector z u . This directly relates the achievable rate to the com¬ 
pressed sensing literature. For a general relation, assume the 
channel fading coherence equals Lc symbols, the effective 
achievable rate of user u (considering the training overhead) 
can be written as 

> Ru (l - j^) (l - (1 - 4 d3) 

where M t equals the number of measurements needed to 
guarantee the support recovery of the sparse vector z u with 
probability at least 1 — e. Given this relation in ( fl3| , it is 
interesting to design the measurement matrix $ to maximize 
this lower bound on the achievable rate. It is also of interest to 
optimize the number of measurements to handle the tradeoff 
shown between the accurate channel estimate (small e) and 
large training overhead. 


5. SIMULATION RESULTS 


In this section, we evaluate the performance of the proposed 
downlink mmWave system operation. We consider the sys¬ 
tem model in Section [2] with the BS employing a ULA of 64 
antennas and 4 MS’s, each with 32 antennas. The system op¬ 
erates at 28 GHz with a bandwidth of 50 MHz. The BS-MS 
distance is 500 m, and all channels are LOS and single-path. 
In the channel estimation phase, the BS and MS’s apply the 
random b eamf orming and measurement procedure described 
in Section 3.1 with Nq S = Nq^ s = 4, and with average trans¬ 
mit power equal to 37 dBm. The beamformers/combiners are 
then built as described in Section [L2l 

In Fig. [2] the per-user achievable rate of the proposed 
system operation is shown versus the number of compressed 


sensing measurements. This figure illustrates that ~ 280 — 
330 measurements are needed to approach the achievable rate 
with perfect channel knowledge. While this number may look 
large, it is actually an order of magnitude less than what is re¬ 
quired by traditional lower frequency solutions, which is ~ 
A r B.s Ams = 2048 symbols. Note also that this training over¬ 
head does not depend on the number of users, which makes 
compressed sensing of special interest to multi-user mmWave 
systems. Compared with adaptive compressed sensing tech¬ 
niques 0 , they may require less training to estimate each user 
channel, e.g., ~ 50 — 100 in |5j. This overhead, however, 
scales with the number of users, which means 200 — 400 for 
4 users, and larger if more users are served. Finally, note that 
these results just illustrate a sufficient number of measure¬ 
ments given the proposed design of the measurement matrix 
$ in Section 3.1 Therefore, optimizing this design may lead 
to even smaller training overhead. 

In Fig. 0 the effective achievable rate is shown for dif¬ 
ferent numbers of fading coherence values, Lc- This figure 
indicates that it is important to compromise between the ac¬ 
curacy of the channel estimate and the training overhead to 
maximize the achievable rate, especially with fast channels. 


6. CONCLUSION 

In this paper, we proposed and evaluated a low-complexity 
downlink system operation for multi-user mmWave systems 
based on compressed sensing channel estimation. Simulation 
results showed that the proposed system operation requires 
a relatively small training overhead (w.r.t. the channel matrix 
dimensions) to achieve a very good performance. Results also 
indicated that the number of measurements need to be wisely 
selected to maximize the effective system sum-rate. 
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